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PAC.E 


KOVEMENT  OF  A WING  WITH  A SFALL  ASPECT  SATIO  NEAR  THE  INTERFACE  OF 
FLUIDS  WITH  DIFFERENT  DENSITIES 


S.  I.  Putilin  (Kiev) 


A.  N.  Panchenkov  [7]  studied  the  froblein  of  the  movement  of  a 
wing  with  a small  aspect  ratio  below  a free  fluid  surface.  The 
integral  equation  for  a wing  with  a small  aspect  ratio  moving  above 
the  interface  of  fluids  with  different  densities  was  obtained  in  a 
similar  manner  in  [H],  But  these  studies  did  net  determine  the 
distribution  of  the  load  over  the  wing  chord.  In  this  report  we  will 
solve  the  problem  which  makes  it  possible  tc  find  the  load 
distribution  over  the  wing  chord  when  condition  XFr*-wO,  is 

satisfied*  the  Froule  number  being  calculated  tor  the  wing  span. 
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1.  The  velcx:ity  potential  ot  an  airfoil  moving  above  a fluid 
interface  should  satisfy  th«^  following  conditions  on  the  interface 
(the  subscript  1 refers  to  the  upper  and  2 - tc  the  lower  fluid): 


v®!,)  - 0; 


where  v is  the  wing  speed.  During  inverse  motion  the  flow 

is  directed  toward  negative  x. 


Using  concept  [5],  by  the  acceleration  potential  method 


j 09  OQ 

K (jc— (^—1))*+ (*  — ;)*"  " ^ J j* 


where  w « (*-_  |)X -p  (y  — tj)*,  , we  obtain  the  following  expression  for  ®r 


dt 


ds. 


(I) 


where 


AT  - —4=.  (V 1/X5+F-. 

2yi.*+k*  ^ ‘ ' 

Q - av  /X*+  k*—  X*+  /liX; 
C*+  Cl 


(2) 

(3) 
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The  roots  of  ejuatioii  Q = 0 lie  in  the  loner  half  plane,  and 
they  must  pass  above  it  duriny  integration  in  the  final  result.  By 
transferring  the  integration  contour  to  the  required  half  plane,  we 
obtain  the  asynptotic  concept  of  the  velocity  potential  at  large 
negative  values  of  •* — 5: 
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The  solution  is 


Yi(l) -6(1-1). 

for  a flat  [-late  viith  a sirall  aspect  ratio  in  a r.  unbounded  fluid.  We 
can  also  anticipate  that  the  majority  of  the  lead  will  be 
concentrated  on  the  leading  edge  tioar  the  interface  of  fluids  with 
ditteren'^  densities  and  we  car.  assume  that 

I 

j"Yi(S)cos[w{Ar  — =2. 

X 


After  this  simplification,  eguaticn  (5)  assumes  the  form  of  the 
equation  for  a wing  with  a finite  span  and  the  eptimum  load 
distribution  solved  by  A.  N.  Panchenkov  and  p.  I,  Zinchuk  [6]. 

Osiny  the  results  obtained  in  this  study,  we  have 


Vtij/)  “ ^ y 1 —!/*(A,+  Ag!/*+Agt/*+  . . .), 
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where  A<  depend  on  the  lead  aiid  are  expressed  ty  spe*cial  functions 
r f^\ 

‘P'f)'  The  expressions  which  detpinine  the  force  are  given  in 

[8  ]. 

2.  A general  expression  for  the  velocity  pcter.*’ial  of  an  airfoil 
moving  above  the  interface  cl  fluids  with  different  densities  was 
obtained  in  [6],  Tuis  expression  can  be  rewritten  as  follows: 
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where  j 


T-  E 

•=  — — cos  0 + (y  — Ti)  sin  6.  (8) 

The  d imeas  ionle:ss  values  found  with  fotaulae  (4)  are  introduced 

here.  Contour  passes  below  particular  point  ^ and 

cos*  6 

contour  Lj  - above  it. 

We  will  study  the  behavior  of  derivative  <Pi  at  A,-*-0.  x<5 
The  first  two  terms  were  studied  in  [2],  Tney  can  be  estiaated  for 
tiiidry  integrals  by  successive  int-'gtation  by  parts  [9],  which 
results  in  this  series: 

- 2 

- J J / (k)  dM  8 « - £ f rm  (z + s ± to)-'-"d9. 

_2  ^ n 

s'  "? 


Here  the  plus  sign  refers  to  contour  Li  and  the  minus  sign  - to 
contour  Lp; 


f(k) 


fe*cos9 
k cos*  0 — (1) 


The  first  non-zero  t^;rm  is  obtained  at  n = 2.  It  is  on  the  order  of 
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cos  Bdd 


2+;±f 


— r-^  COS  8 + — Ti)  sin  8 


Thf  last  term  can  be  festimated  vy  the  stationary  phase  method  [3], 
[9],  For  the  integral 
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the  points  where  H’ (0)  =0  are  aetermiring.  In  cur  case 


//(9) 


1 X(y  — t])  sin  9 
cos  8 X — I cos*8 


and  the  stationary  points  are  determined  by  the  formula 


sin*e=  -1  — Vl  —8d* 
2(1  ’ 
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y— t) 
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where 
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fio.  thp  four  values  of  sinO  U is  necessary  to  select  those 
whose  siqn  is  the  opposite  cf  that  of  cl.  if  we  consider  i to  be  small 
and  disrejard  the  terms  on  the  order  of  cr  hicjher,  we  will  have 


sin’01*  1 _ 4ds; 


sine,=  — i/. 


The  term  of  the  asymptotic  expansion  corresponding  to  root  6,  will 

contain  the  factor  exp  — ] 

Ad*  J 


and  can  be  omitted.  Hoot  corresponds  to  the  term 


(9) 


After  eliminatinj  the  tern-s  whose  crder  of  magnitude 
than  that  of  x from  we  will  have  the  equation 

no, 
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Integration  by  J extends  ovrt  segment 

•P***®  at  with  this  degree  of  accuracy. 

We  will  limit  ourselves  to  the  case  when  the  angle 
does  not  change  over  the  span.  We  will  find  function  Y(£. 
the  toim  of  the  product 

Y (5. 11)  - Y,(g)  Y,(n)- 

Only  cue  of  the  three  terms  in  eguation  (10)  depends  on 
results  in  the  equation 

+1 

J [y  — n (i/— j ” const . 


since 

of  attack 
’l)  in 


y,  which 


(11) 
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Expression  (11)  is  in  the  torin  of  the  equation  for  a winy  with  a 
finite  span  and  the  optimum  load  distribution  moving  above  a screen. 
The  solution  to  this  equation  was  obtained  in  [6]  in  the  form  of  the 
series 

Y,  {y)  = A^+  . .). 

where  the  function  of  At  is  expressed  by  the  parameter 

1 — 2h,  Whence,  setting  the  constant  in  equation  (11)  equal 

to  1 , we  will  have 

+1 


with  consideration  of  this,  equation  (10)  can  be  rewritten  as 

I 1 

^ W - j Y»(5)  + y (1  - a)5*S)  (A)e-^  J y,(1)  x 
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Tho  substitution  of  variables  t),  x-A,(1--/) 

to  the  form 


r / 

J Yi(t)  dx-fi  r Y,(T)  /C  (/  - 1)  dx  - a (/), 


where 


D»_i/C*(p)Y*W-a*(p). 


This  equation  can  bo  solved  by  the  Laplace 
[4].  Indicatinq  the  expression  for  the  function 


have 


Y*W- 


pa*(p)  _ 


reduces  this  equation 

(12) 


(13) 

Carson  transformation 
by  asterisks,  we  will 
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Ihe  oveiall  winy  ch  ai  a c ter  iat  ics  can  be  expressed  directly  by  ; 

thus,  the  problem  will  Le  solved  by  finding  the  resolvent  of  v (t) . 


1.  Similarly,  we  can  study  a winy  itoviny  below  the  interface  of 
two  fluids.  Dsiny  the  results  of  study  [4],  we  will  obtain  the 
following  expression  roc  the  velocity  potential,  analogous  to 
ex  pre  ssi on  (7 ) : 
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The  ccaparisori  of  fotmiiide  (7)  and  (16)  shows  that  the  final 
expression  for  this  case  only  differs  from  e^juation  (12)  by  the  value 
of  which  must  be  replaced  by 

1| 

with  accuracy  up  to  the  terms  on  the  order  of  )J>/3  asymptotic  ■ 

expression  (13)  of  root  K(t)  can  ne  replaced  by  the  expression 

I 

! where  ^ is  the  3essel  function.  Then  we  will  have 

1 
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P,  - 2nfip*—  p*«o  (1  + n*ji*)  + 2nn©V  — 5*; . 

nn/7* — p\i)+rt(io»*p  — (I)*; 

P,- Jifip  (p*+ 


we  win  ha  ve 


v(t) 


-/o(“0+  Pnj 


where  Pi  are  the  roots  of  equation  Pi(p)  =•  0- 

For  an  unbounded  fluid,  = 0 and  v(t)  = 1,  and  we  will  have 

fD(0«a(0. 

If  the  function  a (t)  is  continuous,  (D(/)  is  also  continuous  and 
limited  at  t > 0,  in  particular,  for  the  trailing  edge  of  the  wing. 

This  provides  for  the  satisfaction  of  the  Zhu kc vskiy-Chapl ygin 
condition  [9].  It  is  easy  to  sec  that  the  equation  for  lack  of 
passage  is  not  disturbed  it  we  place  a vortex  of  finite  intensity  T. 
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If  we  set  r— (D(1)y,(t,),  we  will  have 
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